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We discuss phenomenological consequences of requiring the cancellation of quadratic divergences 
up to the leading two- loop order within the Two-Higgs-Doublet Model (2HDM). Taking into account 
existing experimental constraints, allowed regions in the parameter space, permitting the cancella- 
tion, are determined. A degeneracy between masses of scalar bosons is observed for tan /3 > 40. The 
possibility for CP violation in the scalar potential is discussed and regions of tan /3 — M H ± with 
substantial amount of CP violation are determined. In order to provide a source for dark matter in 
a minimal manner, a scalar gauge singlet is introduced and discussed. The model allows to amelio- 
rate the little hierarchy problem by lifting the minimal scalar Higgs boson mass and by suppressing 
the quadratic corrections to scalar masses. The cutoff originating from the naturality arguments is 
therefore lifted from ~ 0.6 TeV in the Standard Model to >2.5 TeV in 2HDM, depending on the 
mass of the lightest scalar. 

PACS numbers: 12.60.Fr, 13.15.+g, 95.30.Cq, 95.35. +d 
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I. INTRODUCTION 



The goal of this work is to extend the Standard Model (SM) such that there would be no quadratic divergences to 
scalar masses up to the leading order at the two-loop level of the perturbation expansion. The quadratic divergences 
were first studied within the SM by Veltman [l[, who showed that applying dimensional reduction Q one gets the 
following quadratically divergent one-loop correction to the Higgs boson (h) mass 



A 2 



(1) 



where A is a UV cutoff and v ~ 246 GeV denotes the vacuum expectation value of the scalar doublet. The issue 
of quadratic divergences was then investigated further adopting other regularization schemes (e.g. point splitting Q) 
and also in Q without reference to any regularization scheme. 

Since precision measurements require a light Higgs boson the correction (|T|) exceeds the mass itself even for small 
values of A, e.g. for m/j = 130 GcV we obtain <$( SM )m| ~ m\ already for A ~ 580 GeV. On the other hand, if we 
assume that the scale of new physics is widely separated from the electro-weak scale, then constraints that emerge 
from analysis of operators of dimension 6 require A> a few TeV. The lesson from this observation is that whatever is 
beyond the SM physics, some amount of fine tuning is necessary; either we tune to lift the cutoff above A ~ 580 GeV, 
or we tune when precision observables measured at LEP are fitted. 1 Tuning both in corrections to the Higgs mass and 
in LEP physics is, of course, also a viable alternative which we are going to explore below. So, we will look for new 
physics in the TeV range which will allow to lift the cutoff implied by quadratic corrections to m\ to the multi-TeV 
range and which will be consistent with all the experimental constraints — both require some amount of tuning. Within 
the SM the requirement 5^ SM ^ ) m 2 l = implies mu ~ 310 GeV. However, as is very well know, the present data favor 
a light SM Higgs boson — according to the PDG after including all the available experimental data and taking 
into account theoretical uncertainties, the 99% CL upper limit for the Higgs mass reads: rrih < 194 GcV. Therefore, 
within the SM the one-loop condition j( SM )m| = requires an unrealistic value of the Higgs boson mass. 

Examining closer the experimental constraints one finds also the following tension which emerges in the pro- 
cess of fitting all the available data to the SM (see @ for a recent review). Hadronic asymmetry measurements 
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(Ap B , Ap B , Qfb) favour a heavy Higgs boson, with rrih ~ 500 GeV, while leptonic asymmetries (Alr, A l FB ) together 
with non-asymmetry precision measurements (mw^Zi • ■ • ) favour a Higgs mass smaller by one order of magnitude. 
If (^-fbj 4fbj Qfb) are omitted from the fit one obtains rrih ~ 50 GeV with an upper limit nrih < 105 GeV at the 
95% CL [gj. Moreover there is the LEP lower limit on the Higgs mass, mh > 114.4 GeV Q- The fit which combines 
all the data is therefore of low quality. That observation suggests a modification of the SM which would allow for a 
heavy Higgs boson with a mass at least above the LEP limit. For that the SM prediction for the oblique parameters 
S and T must be modified by the extension of the SM that we are seeking. 

Here we are going to construct a model which would both soften the little hierarchy problem by suppressing S^ SM 'm h 
and which would allow to lift the central value for the Higgs mass up to a value which is well above the LEP limit 
(presumably it would imply a better fit of the precision observables) . We would like to point out also that increasing 
the Higgs boson mass would ameliorate the little hierarchy problem even if <j( SM )m| was not suppressed (since then 
larger cutoff would lead to the correction of the order of the mass itself). 

There are only two ways to suppress 5^ M ^m\/m\: one can either modify the SM such that (i) larger SM-like Higgs 
boson mass is allowed; or (ii) extra radiative corrections to S^^m^ emerge that partially cancel ([T]). The best-studied 
example of the second approach is provided by supersymmetric theories for which 8m\ <C ra\ up to the GUT scale, 
however the suppression of <j( SM )m| could also be achieved through very modest means, e.g. by introducing just 
extra real scalar singlets to the SM (although more tuning than in the supersymmetric case is necessary). The 
first strategy was followed in Q within the so-called inert doublet model 2 (IDM). There, a second Higgs doublet was 
introduced and an exact 7Li symmetry was imposed to provide a dark-matter candidate. As shown in @ a large 
(400 — 600 GeV) SM-like Higgs boson mass was allowed by the addition of an extra Higgs doublet (the inert doublet). 
It was demonstrated that the extra contributions to the oblique parameters originating from the inert doublet (with 
physical fields H ± , A and S) can cancel large effects of a heavy SM Higgs, such that rrih ~ 400 — 600 GeV is allowed. 
Here we propose a model which does both, i.e. suppression of (5*- SM ^m 2 l by contributions from some extra states (which 
implies reduced i5^ SM ^m 2 ! /m 2 l ) and which modifies the results of the global fit such that much heavier SM Higgs boson 
is allowed (that also helps to decrease 8^ SM ^rn\/rn\ and in addition it could eliminate the tension caused by the high 
LEP lower bound in the presence of a low central value from precision tests). A different approach to this problem 
has been proposed in [l2T | . 

Another well-known problem of the SM is the strength of CP violation (CPV) which is too weak to make the 
electroweak baryogenesis viable jl3j |. It is also worth noting a too slow phase transition (within the SM) which is 
another difficulty for realistic baryogenesis [l3| . 

In light of the above remarks it seems very natural to consider simple extensions of the SM scalar sector, as 
for instance the Two-Higgs-Doublet Model (2HDM). Our intention is to bring the reader's attention to a region of 
parameter space that not only is consistent with standard theoretical requirements (positivity and unitarity) and 
satisfies all the relevant experimental constraints, but also offers a simple pragmatic option to reduce the size of the 
quadratic corrections to scalar two-point Green's functions (so in other words to scalar masses). It has been noticed 
a long time ago [l4|] that within the 2HDM one can cancel quadratically divergent corrections to two-point Green's 
functions for scalar particles. Some phenomenological consequences of the cancellation were discussed already in [l5l ]. 
It is well known [l6[ that within 2HDM extensions the oblique parameters S, T and U can be modified such that SM 
contributions growing with rrih (oc hi rrih) could be canceled by other terms (originating from extra scalars present in 
the 2HDM), so that the lightest Higgs boson could be relatively heavy, see [9[. Within the 2HDM the electroweak 
phase transition could also be made fast enough [l7| to make electroweak baryogenesis viable. The 2HDM provides 
also new sources of CP violation in interactions of neutral scalars. Therefore, here we will discuss 2HDMs which do 
not suffer from quadratic divergences in scalar two-point Green's functions, the tempered Two-Higgs-Doublet Model, 
seeking a model which also allows for CP violation in the scalar potential. Since we have argued above that the heavy 
SM-like Higgs boson would be more consistent with experimental data, we will investigate how much CP violation in 
2HDMs is allowed after lifting the Higgs mass well above the LEP lower limit and by imposing the conditions needed 
to cancel quadratic divergences (so as to ameliorate the little hierarchy problem). We will not address here the issue 
of the electroweak phase transition. 

In a recent publication [l8[, motivated by similar arguments, we have considered a version of the IDM with CP 
violation introduced by replacing the SM-like Higgs doublet by a pair of doublets. There, a candidate for dark matter 
(DM) was provided by the lightest neutral component of the inert doublet (as in the original IDM). In the model 
considered here, CP violation again originates from the 2HDM, however in order to accommodate a DM candidate 
in a minimal manner (instead of introducing the inert doublet as in [l8|) we extend the model by a real singlet. 3 In 
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fact, it is intriguing to note that the singlet is even more inert than the original inert doublet since it interacts only 
with the Higgs doublets and with right-handed neutrinos, having no gauge interactions. 

The paper is organized as follows. In Sec. (jnj we investigate theoretical and phcnomenological consequences of the 
cancellation conditions within the general 2HDM. In order to accommodate a DM candidate we introduce an extra 
real scalar gauge singlet, that is discussed in Sec. IIII1 Section HVl contains our summary. 



II. NON-INERT TWO-HIGGS-DOUBLET MODEL 



A very appealing possibility would be to combine the IDM with the idea of canceling the one-loop quadratic 
divergences. However, as we have shown in poj . that is impossible because ot the vacuum stability conditions in the 
IDM are inconsistent with the requirement of cancellation of quadratic divergences. Since our intention is to build a 
2HDM, which has no quadratic divergences at least at the one-loop level therefore, in the following we will consider 
a general (non-inert) 2HDM hoping for both a successful implementation of the cancellation condition and for a new 
source of CP violation. The price to pay will be the loss of a DM candidate. We return to that issue in Sec. IIIII 

In order to accommodate CP violation we consider here a non-inert 2HDM with softly broken 7Li symmetry which 
acts as $1 — > —$1 and ur — > —iir (all other fields are neutral). The scalar potential then reads 



V ((/>i, </> 2 ) = {w?!^! + 77122^ 



22^2^2 



l 12Vl l 



H.c. 



~Ai(^i) a + ^a 2 (40 2 ) 2 + A 3 (</) t i0 1 )(40 2 ) + a 4 (4^)(4^: 



(2) 



A 5 (4<M 2 + H.c. 



The minimization conditions at 



6?) = Vl /y/2 and ($} = v 2 /y/2 can be formulated as follows: 

= vf\i + u|(A 3 4 5 - 2i/), 
= W2A2 + w 2 (A 3 45 - 2u), 
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'22 



(3) 



where A345 = A3 + A4 + Re A5 and v = Rem 2 2 / (2vii) 2 ). 

We assume that 4>i an d 4>2 couple to down- and up-type quarks, respectively (the so-called 2HDM II). 



A. One-loop quadratic divergences 

The cancellation of one-loop quadratic divergences for the scalar two-point Green's functions at zero external 
momenta {G%, i = 1, 2) implies [14[ in the case of 2HDM type II: 



3 2 


3 2 


v 2 




— m 7 - 
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v 2 


= 2 m w ~ 




h Y 



3, . 1 
2 Al + A3+ 2- 

3, . 1 



2- - 0, (4) 

G 2 = + - A m 2 z + - ^A 2 + A 3 + - A 4 J - 3^f = 0, (5) 

where v 2 = vf + v 2 , tan j3 = w 2 /wi and we adopt the notation: sp = sin /3 and cp = cos j3. We note that when tan/3 is 
large, the two quark contributions can be comparable. In the type II model the mixed, <p\ — 4>2, Green's function is 
not quadratically divergent. 

In the general CP-violating case, the quartic couplings A^ can be expressed in terms of the mass parameters and 
elements of the rotation matrix needed for diagonalization of the scalar masses (see, for example, Eqs. (3.1)-(3.5) of 
EH): 

Ai = -^[c\c 2 2 Ml + (c lS 2S 3 + sic 3 ) 2 M 2 2 

+ (c 1 s 2 c 3 ~s 1 s 3 ) 2 M 2 ~sj f i 2 }, (6) 
2^[s 2 C^M 2 + (C1C3 - Sl s 2 S3) 2 M 2 2 
2 .,2 



S/3 v 



+ (c 1 s 3 + s 1 S2C 3 ) 2 M 2 -cj ! n}, (7) 
{c lS i[c 2 M 2 + ( S 2 S 2 -c 2 )M 2 



CpSpV 2 
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+ (44 - s 2 )M 2 ] + s 2 c 3 s 3 (c 2 - s\){Ml - M 2 )} + i[2M| ± - M 2 ], (8) 

A 4 = 1 [ S 2 M 2 + C 2 S 2 M 2 + C 2 C 2 M 2 + ^2_ ^2^ ] ^ (g) 



Re A 5 = ^ ~ c\s\Ml - 44Mi + /i 2 ] , (10) 

ImA 5 = — x{cp[cic 2 s 2 M 2 - c 2 s 3 (as 2 s 3 + sic 3 )M 2 2 

+ c 2 c 3 (sis 3 - cis 2 c 3 )Ml\ + s / j[sic 2 s 2 M 1 2 (11) 

+ C 2 S3(ciC3-S 1 S2S3)M 2 2 -C2C 3 (ciS3 + SlS2C3)M 3 2 ]}, 

where /i = v 2 v while Cj = coscti and Sj = sin a; refer to the neutral- Higgs-sector rotation matrix R, the latter 
parametrized in terms of the angles a\ , a 2 and a 3 according to the convention of [22[ . 

It will be useful to adopt the following relation (emerging from the diagonalization of the neutral Higgs mass 
matrix [Hj]) between M 2 , Mf and Mf: 

M 2 = M 2 R 13 (-R n + fl 12 tan ,3) + MjR 23 {~R 21 + i? 22 tan/3) 
3 R 33 (R 31 - R 32 tan (3) ' 1 j 

Substitution of (|T2"j) into ([6|)- (fTT|) allows to express the quartic couplings through the mixing angles together with 
M 2 , Mf , Mjt ± and /x 2 (eliminating M|). Then, inserting the appropriate quartic couplings into the conditions for 
cancellation of the quadratic divergences (J4J) d5j) we obtain two linear equations for Mf and Mf with coefficients 
depending on the mixing angles on (as well as on M^ I± and fi 2 ). Therefore, for a given choice of ctj's, the squared 
neutral-Higgs masses Mf, M 2 and Mf can be determined from the cancellation conditions (HJ)-© in terms of tan /3, 
fj, 2 and Mfj ± . 

Scalar masses resulting from a scan over c^, Mjj± and tan (3 are shown in Fig.[T]for /j, = 200 GcV and // = 500 GcV. 
The charged Higgs boson mass was varied between 300 GeV and 700 GeV. Since only large tan/3 will turn out to 
be allowed we have chosen to display plots with 40 < tan/3 < 50 in order to illustrate a specific property of the 
scalar spectrum that is visible at large tan (3. Under the scan, the M 2 were calculated along the lines described 
above. The only extra constraints (the cancellation of quadratic divergences was, of course, guaranteed implicitly by 
the construction) imposed were Mf > and Mi < M 2 < M 3 . A striking degeneracy of the neutral-Higgs masses 
is observed for the case of large tan/3. This degeneracy can be understood by expanding M 2 for large tan/3. The 
cancellation conditions, Eqs. (0])-(© can then be expressed as follows: 

Y n M 2 + Y 12 M 2 - Y 13 (4m 2 + A ; 2 ) = O ( -^—] , (13) 



tan /3 
1 

tan j3 



where 4 



2i? 12 i? 22 tan/3(Mf - Mf ) - i? 33 [-4m 2 - 2M 2 ± + 12m 2 + f i 2 } = O ( — ^ ) , (14) 



Yii = — R\ 2 R\ 3 R\i + R11R32R33, 

?2 1 r>2 



Y12 — — R22R23R31 + -^21-^32^33, 

Y13 = -R32-R331 

— 2 ^ 2 ^ 2 / -i r- \ 

to = -m w + -m z . (15) 
First, it is useful to notice that (fT4"|) implies 

M 2 - M 2 2 ~ - 1 R ™ [-Am 2 - 2M| ± + 12to 2 + (16) 
tan/3 R\ 2 R 22 

Therefore, for tan/3 3> 1, we expect to have Mf ~ Mf unless |i?i 2 i? 22 | < 1 or — 2M%± + /x 2 is very large. 5 Secondly, 
solving ([T3| -(fr4 |) one finds that to leading order (for large tan 0) Mf = Ml — fi 2 + Am 2 . 



4 Note that Y u + Y 12 = Y 13 . 

5 Note that cancellations between the and /x 2 terms are possible. This is why the degeneracy survives even for fj, as large as 
fx = 500 GeV, see Fig. [TJ 
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FIG. 1: Distributions of allowed masses M% vs Mi (left panels) and M3 vs M2 (right), resulting from a scan over the full range 
of an, tan/3 G (40,50) and M H ± £ (300,700) GeV, for fi = 500 GeV (top) and ^ = 200 GeV (bottom). No constraints are 
imposed other than the cancellation of quadratic divergences (|4}-([5}, Mf > and Mi < M2 < M3. The color coding indicates 
increasing density (while scanning over the parameter space) of allowed points as one moves inward from the boundary. 



On the other hand, expanding (jT2|) for tan/3 ^> 1 one obtains 



r 2 _ M 1 Ri2Ri3 + M2R22R23 

R32R33 \tan/3 



Therefore (invoking unitarity of i?) it is seen that the degeneracy Mi = M2 implies that also Mi = AI2 = M3. 
Finally we can conclude that for large tan/3 one obtains Mi ~ M 2 ~ M3 ~ // 2 + 4?tz|, this explains the approximate 
degeneracy observed in Fig. [T] 6 



The reader should be warned that the above expansions are justified if the coefficients of sub-leading terms cx 1/ tan /3 are not enhanced 
by special values of the mixing angles (that would correspond to CP conservation in the scalar sector). Since here we are interested in 
the case of CP violation, we will not elaborate on those CP conserving limits. 
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B. Two- loop leading quadratic divergences 

The generic form of the quadratically divergent contributions to scalar two-point Green's functions at zero external 
momenta reads Q 



SG, = A 2 £/W(A) 



In I ^ 



A x 1 " 



(18) 



where n corresponds to (n + l)-loop contribution, A stands for relevant coupling constants, fi is the rcnormalization 

scale and fn\x) is a calculable (order by order) function (polynomial) of the couplings. It should be noticed that at 
the (n + l)-loop level there exist also sub-leading contributions that contain terms oc A 2 (ln A) m with m < n, so for 
instance at the two-loop level the leading contribution is oc A 2 In A while there are also sub- leading terms oc A 2 . The 
coefficients of the leading terms, /«, (A), can be determined recursively adopting a nice algorithm noticed by Einhorn 
and Jones 0]: 

(n + l)/# a = t§=ti ] =E^^7^ ( 19 ) 

where the sum runs over coupling constants that contribute to the coefficient fn- Hereafter we will limit ourselves 
to the leading two-loop contributions. Therefore, to calculate f[ l \ only the one-loop coefficient and one-loop 



beta functions are needed. As beta functions for the 2HDM arc known [24[ the cancellation condition for quadratic 
divergences up to the leading two-loop order can easily be determined: 



d + Sd = and G 2 + 5G 2 = (20) 



with 



8G l = ^ [9 52j S 92 + 3< 7l /3 Sl +6/3 Al +4/3 A3 +2/3 A4 ] In (j) (21) 

5G 2 = V - [9g 2 [3 g . 2 + ZgiP 9l + 6/?a 2 + 4/? A3 + 2/3 A4 - 24g t f3 9t ] In (22) 

In what follows, adopting (|B|)- ([TTT) we will be solving the conditions (j2U)) for the scalar masses Mf for a given set of 
ai's, tan/3, [i 2 and M^ ± . For the renormalization scale we will adopt v, so ft = v. 7 Then those masses together with 
the corresponding coupling constants, will be adopted to find predictions of the model for various obscrvablcs which 
can be confronted with experiments. 

C. Positivity and unitarity constraints 

The requirements of positivity for the 2HDM model potential are well known (Io| : 

Ai, 2 > 0, (23) 
A 3 >-v / AiA 2 ", Xl = A3 + A4- |A 5 | >- V / AiA 2 ". (24) 

One could also require the above conditions to be satisfied up the unification scale A ~ 10 15 GeV. That approach, 



resulting in much stronger constraints in terms of allowed scalar masses and tan /3 was followed in Ref . jll| . Here, as 
we consider UV completion appearing at the scale of a few TeV we do not follow that line of reasoning. 



Since we are using tree-level relations between quartic couplings and scalar masses, the renormalization scale should be of the order of 
the masses themselves, that is why we adopt here p, = v. For a more exhaustive discussion of the renormalization-scale dependence, see 
the first paper of [25| . 
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FIG. 2: Two-loop allowed regions in the tan/3-M ff ± plane, for A = 2.5 TeV, for (i = 300,400,500 GeV (as indicated). Red: 
positivity is satisfied; yellow: positivity and unitarity both satisfied; green: also experimental constraints satisfied at the 95% 
C.L., as specified in the text. 
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FIG. 3: Similar to Fig. [2] for A = 6.5 TeV. 



D. Experimental constraints 

We impose the following experimental constraints: 

• The oblique parameters T and S 

• Bq — Bq mixing 

• B -> X sl 

• B -> tv t X 

• B -> _Dr^ r 

• LEP2 Higgs-boson non-discovery 

• Rb 

• The muon anomalous magnetic moment 

• Electron electric dipole moment 

For more details concerning the implementation of the experimental constraints, see refs. [3 [U |3(|. Subject to 
all these constraints, we find allowed solutions of (J^QJ) - The recent paper (2(| also contains an exhaustive analysis 
of experimental constraints on the 2HDM type II. The lower limit on the charged Higgs-boson mass adopted here, 
M T > 300 GeV (basically determined by the b — > 57 constraint) agrees roughly with the 95% CL lower limit, 316 GeV, 
obtained in [2(| irrespectively of tan j3. 
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E. Allowed regions 



Imposing the above conditions we find allowed regions in the tan /3-M#± plane as illustrated by the red domains 
in the tan/3-M#± plane, see Figs. [2] and [3] for fixed values of \x. The allowed regions were obtained scanning over 
the mixing angles a.j and solving the two- loop cancellation conditions (|20[) . Imposing also unitarity in the Higgs- 
Higgs-scattering sector [27l - l29| (yellow regions), the allowed regions are only slightly reduced. Requiring that also 
experimental constraints listed in the Sec. HID I are satisfied one obtains green regions shown in the figure. 

For parameters that are consistent with unitarity, positivity, experimental constraints and the two-loop cancellation 
conditions (|2T)]) . we show in Figs. HJ[S] scalar masses resulting from a scan over a 2 ;, M H ± and tan/3. Those plots could 
be compared with Fig. [T] One should however remember that in the two-loop case also unitarity, positivity and 
experimental constraints are taken into account. Note that in Figs. [2] and [3l consistent solutions are obtained only 
for tan/3> 15 while at the one-loop level, also a small low-tan /3 region was allowed after imposing all the constraints. 
That small low-tan /3 region is disallowed after the two- loop corrections are imposed, see [20j for the one- loop result. 
As we have noticed for the one-loop spectrum, large tan/3 implies similar scalar masses. This is indeed what is being 
observed in Figs. 2][5] also for the two-loop case. The allowed solutions "peak" around M H ± ~ /i with 20 < tan/3< 50. 
For fj, = 200 and 600 GeV there are hardly any solutions for A = 2.5 TeV and no solutions were found for A = 6.5 TeV. 



F. CP violation 



Here we are going to discuss the possibility of having CP violation in the scalar potential @, subject to the two- 
loop cancellation of quadratic divergences (|20[) . In order to parametrize the magnitude of CP violation we adopt 
the rephasing invariants introduced by Lavoura and Silva [3l[ (see also (32J). We shall here use the basis- invariant 
formulation of these invariants J\ , J 2 and J3 as proposed by Gunion and Haber [33j . As is proven there (theorem 
#4) the Higgs sector is CP-conserving if and only if all Ji are real. In the basis adopted here the invariants read [lH : 



2„,2 



ImJj = -^(Ai -A 2 )ImA 5 , (25) 

V 

2 2 

ImJ 2 = [((A 1 -A 3 -A4) 2 -|A 5 | 2 )« 1 +2(A 1 -A 2 )ReA5^« 2 2 

- ((A 2 - A 3 - A 4 ) 2 - |A 5 | 2 ) v 4 2 ] ImA 5 , (26) 

2 2 

ImJ 3 = ^ 2 -(A 1 -A 2 )(A 1 + A 2 + 2A 4 )ImA 5 . (27) 

It is seen that there is no CP violation when LrnAs = 0, see [l8[ for more details. 

As we have noted earlier, tan/3 above ~ 40 implies approximate degeneracy of scalar masses. That could be 
catastrophic for CP violation since it is well known that the exact degeneracy Mi = M 2 = M3 results in vanishing 
invariants Im Ji and no CP violation (exact degeneracy implies ImAs = 0). Using the one-loop conditions (d))-© one 
immediately finds that A x — A 2 = 4(m 2 /c 2 j — m 2 /si)/v 2 , which implies 



ImJ 1 ^4ImA 5 p " 2 p =-4ImA 5 ^ lO (28) 
rr V v J \ tan / 

In fact the above result shows even more than we have anticipated. If tan j3 is large then Im Ji is suppressed not only 
by ImAs ~ (as caused by Mi ~ M 2 ~ A/3) but also by the factor (ml/v 2 ), as implied by the cancellation conditions 
@-([S])- The same suppression factor appears for Im J3. The case of Im J 2 is more involved, however when m 2 /v 2 is 
neglected all the invariants (|25l (|27|) have the same simple asymptotic behavior 

for large tan/3. It is also worth noticing that tan/3 = m t /mb(— 38) implies Ai = A 2 , which in turn leads to exact 
vanishing of Im J\ and Im J3. Qualitatively those conclusions survive at the two- loop level. For a quantitative 
illustration we plot in Figs. |B][7] maximal values of the invariants in the tan/3— M#± plane with all the necessary 
constraints imposed, seeking regions which still allow for substantial CP violation. At high values of tan/3 these 
invariants are of the order of 10 -3 , in qualitative agreement with the discussion above. Note that the SM corresponding 
invariant ImQ = {V u dV c bV* b V* d ) ~ 2 x 10~ 5 siuSkm [13, for Vij and 5km being elements of the CKM matrix and CP- 
violating phase, respectively Therefore the model considered here offers at least two orders of magnitude enhancement 
comparing to the SM. 
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FIG. 4: Two-loop distributions of allowed masses M2 vs Mi (left panels) and M3 vs M2 (right) for A = 2.5 TeV, resulting 
from a scan over the full range of on, tan/? G (0.5, 50) and M H ± £ (300, 700) GeV, for (i = 300,400, 500 GeV. Red: Positivity 
is satisfied; yellow: positivity and unitarity both satisfied; green: also experimental constraints satisfied at the 95% C.L., as 
specified in the text. 



G. Stability and the determination of the cutoff 

It should be emphasized here that the conditions (|20l) eliminate the quadratic divergences only up to the leading 
two-loop corrections. Even though the sub-leading two-loop and higher effects are suppressed by powers of coupling 
constants and powers of l/(167r 2 ), nevertheless since the In A term is growing, there exists always A large enough, 
that the hierarchy problem reappears: loop corrections to masses are again of the order of the masses itself. In fact 
that observation allows to determine the value of the cutoff up to which higher order corrections do not reintroduce 
the hierarchy problem (see [1^] for the analogous strategy within the SM) . In general, quadratic corrections to scalar 
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FIG. 5: Similar as Fig. g] for A = 6.5 TeV for fi = 300,400, 500 GeV. 



masses have the form of (fTSl ) 



SMf = A 2 £/£>(A) 



v 



where v is chosen as a renormalization scale. The following naive estimation of /„ is sufficient 
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where the relevant coupling constants were conservatively assumed to be of the order of An. 8 Here we choose as the 
cutoff the maximal value of A such that the higher order corrections do not exceed the mass of the lightest scalar, 



A<4ttMi 

200 (500) GeV the cutoff is at least at A 



(32) 

2.5 (6.3) TeV. Of course, larger Mi would imply 



Then, e.g. for M\ 
higher A. 

Having the cutoff determined, we should address the issue of higher-loop corrections to the equations (|20| that 
ensure vanishing of the quadratic corrections up to leading two-loop effects. As is seen from ([3TJ|) . the generic form of 
the condition for vanishing quadratic divergence is the following 



A 



A- 



(4tt) 2 (4tt)' 



In I* 



A 2 



A 3 



v J (4tt) 4 (4?r) e 



= o. 



(33) 



where A stands for a typical coupling constant. The last two terms shown above (sub- leading two- and leading three- 
loop effects) have been neglected in the present analysis. It is then easy to see that even for the cutoff as large as 
A = 6.5 TeV using a very conservative (large) value for the typical coupling, A = 47r, the precision of the adopted 
approximation is of the order of 12%. Note that whenever A < 47r or A < 6.5 TeV, the adopted approximations work 
better. 



III. 2 DOUBLET + 1 SINGLET HIGGS MODEL: THE CASE FOR DARK MATTER 

In this scenario we combine CP violation present in the non-inert 2HDM (allowing for softly broken W12 symmetry) 
with a real scalar 92 which is a gauge singlet. The singlet provides a natural DM candidate (see [HJ], [35| and Q). In 



This estimate agrees qualitatively with the two-loop result obtained for /1 in the SM, see Eq. (21) in [25ll . 
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this case the scalar potential is the following 

V ((/>i, (/> 2 ) = --\m\ r <j)\(j) 1 +ml 2 4)\(j)2 + 



m\ 2 (f>\(j)2 



H.c. 



•>2<P2) 



+A 3 (<^i)(4<fe) + A 4 (0 t 1 2 )(^0 1 ) + - [\ 5 {<f>\ct> 



H.c. 



(34) 



Note that the term oc p 2 (j>\cf>2 is forbidden as it breaks the 7L 2 symmetry in a hard way. Since tp is supposed to be 
the DM candidate, in order to ensure its stability we have imposed an extra discrete symmetry 7L' 2 such that tp — > —ip 
while other fields are neutral. The symmetry excludes terms odd in tp. The potential should be arranged such that 
the symmetry remains unbroken, so that (tp) = 0. For that it is sufficient to require 



nl > & A„, 771, 772 > 



(35) 



Then it is easy to see that if the standard 2HDM stability conditions (f23|) -([24 |l are fulfilled then the potential ([34| is 
also positive definite. For the mass of the singlet we obtain: vr? = 2fx 2 + r\\v\ + r\iv\. 

Since tp is a gauge singlet, therefore in the presence of right-handed neutrinos which are also gauge singlets the 
following Yukawa interaction is allowed @ : 



£y = -tp(vn) c Y v v R + H.c. 



(36) 



Note that for a number of right-handed neutrino flavours greater than 1, the Yukawa matrix Y v is in general (depending 
on the quantum numbers of vr under Z 2 , see Q) non- vanishing. 

In this model the conditions for cancellation of quadratic divergences are slightly modified: 



i 2" 1 



^A 1 + A 3 + iA 4 ) -3^ = 0, 



u 2 /l 3, x 1 
tU 7?2+ 2 A2 + A3+ 2 A4 
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%+4(r 7l + r ?2 )-8Tr{y ¥ ,Ft } 



0. 



(37) 



The last condition above guarantees vanishing quadratic divergence in corrections to the tp mass. Since for the 
positivity we assumed \ v ,r)i,r)2 > 0, it is clear from the above equation that the presence of the Yukawa coupling Y v 
is mandatory to extend the condition for cancellation of quadratic divergences to the singlet field as well. It should 
also be mentioned that the presence of the singlet docs influence the two-loop corrections to the quadratic divergences, 
those effects are neglected as being small (oc ?/,). 

As we have already mentioned the extra singlet ip provides a candidate for the DM. To estimate its present abundance 
we consider the dominant annihilation channels for (p. The Lagrangian describing relevant cubic and quartic scalar 
interactions reads 



C = -ip\n % vH % + A,,//,//, + \ ± H + H~) 



(38) 



where 



niRiiCp + r] 2 Ri 2 S/3 

1 

2 



CgRisRja) 



Ah 



Vis} 



(39) 
(40) 
(41) 



A detailed study of the DM within this model (with extended conditions for the cancellation of quadratic divergences 
(1371) ) will be presented elsewhere [36| , However here we would like to show that it is indeed natural to expect the 
right DM abundance in the presence of the singlet. For an illustration we will assume that the DM annihilation cross 
section is of the order of the contributions from the lightest neutral Higgs boson Hi of mass M\. 

For an estimate of the DM abundance, we will consider two (pip annihilation mechanisms. First we assume that ipip 
annihilate to 77, qq, , W + W~ and ZZ through s-channel Hi exchange. Then, following [38[ we obtain 9 in the 
non-relativistic approximation the following result for the thermally averaged annihilation cross section 



2„,2 



An\v 



{Am 2 , 



Ml) 2 



T Hl (2m v 



2m u 



(42) 



where Th 1 (2m v ) stands for the decay width of Hi calculated at Mi = 2m v (in the following numerical calculations we 
will use the SM width for the estimate). Now we have to add the contribution from the H x Hi final state. There are 
two contributions: due to s-channcl Higgs exchange, and due to the four-point coupling. We find in the non-relativistic 
approximation 



1 1 

32-7T m 2 , 



Ml 

2 



1/2 



6(m v - Mi) 



An 



KlAmW 2 



Am 2 ~ M 2 + iMiT Hl 



(43) 



where the quartic (p 2 HiHi coupling An is defined by Eq. (|40j) and the trilinear H1H1H1 coupling normalized to v 
is denoted by Am :39j. For an order-of-magnitude estimate of the DM abundance, we will use here Am = 3Ml/v 2 
(this choice, together with (|44p. reproduces results which would be obtained for the SM Higgs doublet tfisM coupled 
to the singlet through the term rjtp \</>sm\ ) an d parameterize Ki and An through one variable 77 as follows: 



2An — Ki 



(44) 



Then, following [381 ] for cold relics one has to solve the following equation to determine the freeze-out temperature 



from Xf 



,/Tf. 



x f 



In 



(45) 



This cross section is also to be found in the literature [35l . however our result is smaller by a factor of 2. We have included both the 
combinatoric factor 1/2 in ipipH\ vertices and statistical factors (both in the initial and in the final states) in (crv). 
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FIG. 8: Inert-scalar coupling r\ (vs m v ) required by the observed DM abundance CIdmIi 2 = 0.106 ± 0.008 [37 
band. As indicated above each panel, the lightest Higgs-boson mass ranges from Mi = 100 to 400 GeV . 



within a 3-<j 



whcrc (av) = (<Jv)i + {uv)2 and g* counts rclativistic degrees of freedom at annihilation and mp\ denotes the Planck 
mass. It turns out that in the range of parameters we are interested in, Xf ~ 0(25), so that this is indeed the case of 
cold dark matter, it also implies that g+ ~ 10 — 100. Then the present density of y's is given by 



n v h 2 = i.07- io 9 — 



x f 



(46) 



g+' m Pi {<Tv) 



In Fig. [8]we show the 3-cr allowed band in log(Ty) vs. m^, as constrained by the observed DM abundance Qrjjvf/i 2 = 
0.106 ± 0.008 (37j . For rrih = 100 and 200 GeV we observe consequences of resonant behavior at m v = The 
thresholds seen at m v ~ 25 and 80 GeV are caused by the rapid change in g+ as a function of temperature and by 
the opening of the W + W~ channel for the decay of a Higgs boson of mass m/i = 2mw, respectively. 

One can conclude that the singlet could indeed provide a realistic candidate for DM: for any m v between ~ 1 GeV 
and ~ 500 GeV there exists an allowed value r\ for which tt v h 2 agrees with the experimental data. Note that if we had 
found only solutions with 77 > 1 and light if (m v <v) then this scenario would be jeopardized since the minimization 
condition requires m 2 > r\\v\ +r]2V 2 (as /i 2 > 0). As seen from Fig. El this is not the case. Most of the allowed region 
corresponds indeed to 77 < 10" 1 if the singlet mass is not too low. 



IV. SUMMARY 



The goal of this work was to build a minimal realistic model which would allow for softening the little hierarchy 
problem through suppression of the quadratic divergences in scalar boson mass corrections and through lifting the 



15 



mass of the lightest Higgs boson. That could be accomplished within Two-Higgs-Doublet Models. Phcnomcnological 
consequences of requiring no quadratic divergences in corrections to scalar masses within the 2HDM were discussed. 
The 2HDM type II was analyzed taking into account existing experimental constraints. Allowed regions in the 
parameter space were determined. An interesting scalar mass degeneracy was observed for tan/3>40. The issue of 
possible CP violation in the scalar potential was addressed and regions of tan j3 — M H ± with substantial strength of 
CP violation were identified. In order to accommodate a possibility for dark matter a scalar gauge singlet was added 
to the model. Requirements necessary for correct present abundance of dark matter were estimated. 

The model we considered here allows to soften the little hierarchy problem by lifting the minimal scalar Higgs boson 
mass and by suppressing the one-loop quadratic corrections to scalar masses. The cutoff implied by the naturality 
arguments is lifted from ~ 600 GeV in the SM up to at least >2.5 TeV, depending on the mass of the lightest scalar. 
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